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Circular-cylinder  model*  having  nominal  diameter*  of  1/8,  1/4,  Circular-cylinder  models  having  nominal  diameters  of  1/8,  1/4, 

and  1  inch  were  mounted  tran*ver*ely  to  the  flow  at  the  nozzle  and  1  inch  were  mounted  transversely  to  the  flow  at  the  nozzle 

centerline.  A  hemisphere  cylinder  having  a  nominal  diameter  of  centerline.  A  hemisphere  cylinder  having  a  nominal  diameter  of 

1  JZ  inch  wa*  mounted  at  the  nozzle  centerline  with  it*  axi*  1/2  inch  was  mounted  at  the  nozzle  centerline  with  its  axis 
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FOREWORD 


The  study  on  which  this  report  is  based  constitutes  a  part  of  the 
research  program  on  aero- the rmo-characteristics  of  flight  of  hypersonic 
vehicles  under  the  U.S.  Air  Force  Contract  AF  33(6 16)-60Z5,  Project  No. 

7064  "Thermo- Aerodynamic  Characteristics  at  Hypersonic  Mach  Numbers," 
Task  No.  70169  "Research  on  Aerodynamic  Flow  Fields.  "  The  project  is 
administered  by  the  Fluid  Dynamics  Research  Laboratory,  Air  Force 
Research  Division,  Wright- Patter  son  Air  Force  Base,  with  Mr.  Robert 
D.  Stewart  acting  as  Task  Scientist. 

As  indicated  in  the  text,  a  portion  of  the  work  (the  hemisphere-cylinder 
experiments  and  the  theoretical  work  of  H.  K.  Cheng)  was  performed  under 
Contract  Nonr  2653(00)  with  the  U.S.  Navy  Office  of  Naval  Research. 

Some  of  the  results  contained  in  this  report  (the  transverse- cylinder 
data)  were  presented  in  a  talk  entitled  "Low-Density  Research  in  the  Hyper¬ 
sonic  Shock  Tunnel"  given  at  the  "Second  International  Symposium  on  Rarefied 
Gas  Dynamics"  held  at  the  University  of  California,  Berkeley,  California 
from  August  3  to  6,  I960.  A  written  version  of  the  talk  will  be  published 
in  the  proceedings  of  the  Symposium. 

The  authors  wish  to  thank  their  colleagues  in  the  Aerodynamic  Research 
Department  for  their  contributions  to  this  program,  particularly  Dr.  H.  K. 
Cheng  for  the  many  consultations  regarding  his  theoretical  analysis.  They 
would  also  like  to  acknowledge  the  assistance  of  Miss  Eileen  Stager  who 
performed  the  detailed  data  reduction  computations. 


ABSTRACT 


A  study  of  stagnation-point  heat  transfer  to  two-  and  three-dimensional 
bodies  at  low  Reynolds  numbers  has  been  undertaken  at  the  Cornell  Aero¬ 
nautical  Laboratory,  Inc.  The  purpose  of  this  investigation  was  to  obtain 
hypersonic  low  density  data  and  to  demonstrate  the  suitability  of  the  hyper¬ 
sonic  shock  tunnel  as  a  facility  for  research  in  rarefied  gasdynamics. 

Low-density  experiments  on  a  transverse  cylinder  and  a  hemisphere 
cylinder  have  been  made  in  air  at  Mach  numbers  from  8.  4  to  1 1.  6  in  the 
CAL  11-  by  15-inch  hypersonic  shock  tunnel.  The  stagnation  pressure  and 
temperature,  as  well  as  model  diameter,  were  varied  to  obtain  a  range  of 
Reynolds  numbers  from  1 1  to  1000,  based  on  flow  conditions  behind  the  bow 
shock  wave  and  model  radius. 

Circular -cylinder  models  having  nominal  diameters  of  1/8,  1/4,  and 
1  inch  were  mounted  transversely  to  the  flow  at  the  nozzle  centerline.  A 
hemisphere  cylinder  having  a  nominal  diameter  of  1/2  inch  was  mounted  at 
the  nozzle  centerline  with  its  axis  aligned  with  the  flow.  Thin-film  resist¬ 
ance  thermometers  were  mounted  along  the  stagnation  line  to  measure  the 
stagnation-point  heat-transfer  rate  at  one  or  more  points  along  the  cylinder 
and  at  the  stagnation  point  on  the  hemisphere.  The  experimental  heat- transfer 
rates  are  presented  and  compared  with  theoretical  predictions.  The  data 
obtained  with  the  transverse  cylinders  are  in  good  agreement  with  continuum 
boundary- layer  theory  at  all  but  the  lowest  Reynolds  numbers  and  highest 
Knudsen  numbers.  The  hemisphere -cylinder  data  indicated  the  presence  of 
a  significant  vorticity-interaction  effect.  Good  agreement  with  the  theoretical 
analysis  of  H,  K.  Cheng  accounting  for  this  effect  is  indicated. 

The  success  of  these  experiments  in  a  shock  tunnel  that  was  not  designed 
for  low-density  operation  has  led  to  a  fuller  investigation  of  the  application  of 
the  shock  tunnel  to  research  studies  of  hypersonic,  rarefied  gasdynamics. 

The  results  of  these  considerations  are  discussed. 
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INTRODUCTION 


The  phenomena  encountered  in  hypersonic,  rarefied-gas  flows  recently 

have  been  the  object  of  an  increasing  number  of  investigations.  Much  of  this 

interest  has  been  stimulated  by  the  prospect  of  manned-hypersonic  vehicles. 

Because  of  the  glancing  re-entry  trajectories  and  the  requirements  for 

sustained  hypersonic  flight,  these  vehicles  will  encounter  rarefied- gasdy- 

namic  phenomena  when  they  are  flying  at  very  high  hypersonic  speed  and 

high  altitudes.  Hypersonic  rarefied-gas  flows  have  been  treated  theoretically 

in  recent  years;  however,  experimental  low-density  data  has  been  confined 

mostly  to  Mach  numbers  of  6  or  less  (see  for  example,  References  1, 2,  3). 

Because  of  its  operating  capability  over  a  wide  pressure  range,  it  was  felt 

that  the  shock  tunnel  could  be  used  to  obtain  hypersonic  low-density  data. 

Consequently,  a  program  of  rarefied-flow  studies  has  been  undertaken  at 

the  Cornell  Aeronautical  Laboratory,  Inc.  utilizing  the  existing  11-  by  15- 

4  5 

inch  hypersonic  shock  tunnel  1  .  The  purpose  of  this  program  was  twofold; 

first,  to  obtain  hypersonic  low-density  data,  and  second,  to  investigate  the  . 
shock  tunnel  as  a  hypersonic  low-density  research  facility. 

Stagnation-point  heat  transfer  to  two-  and  three-dimensional  bodies 

was  selected  as  the  subject  of  the  investigation  both  because  it  is  of  practical 

7 

interest  and  has  been  treated  theoretically  by  Probstein  and  Kemp  ,  Herring  , 

8  9  10 

Hoshizaki  ,  Ho  and  Probstein7  and  Ferri,  Zakkay  and  Ting  .  These  theories 

indicate  that  the  stagnation-point  heat-transfer  rate  will  be  greater  than  that 

predicted  by  boundary-layer  theory. 

Transverse  cylinders  having  nominal  diameters  of  1/8,  1/4  and  1  inch, 
and  a  hemisphere  cylinder  having  a  diameter  of  1/2  inch  have  been  tested. 

*  11 
The  transverse- cylinder  data  were  initially  reported  by  the  authors  at  the 

Second  International  Symposium  on  Rarefied  Gas  Dynamics  held  at  the 

University  of  California,  Berkeley,  California  from  August  3  to  6,  I960, 

and  will  be  published  in  the  proceedings  of  the  Symposium. 


Manuscript  released  by  the  authors  December  15,  I960  for  publication  as 
an  ARL  technical  report.  ' 
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The  experiments  were  made  in  air  at  Mach  numbers  from  8.  4  to  11.6  and 

Reynolds  numbers  from  1 1  to  1000  based  on  flow  conditions  behind  the  bow 

shock  wave  and  model  radius.  The  stagnation  temperatures  ranged  from 

1500  to  2200°K  and  from  3200  to  3700°K.  The  Knudsen  number  based  on 

-4 

mean  free  path  behind  the  shock  and  body  radius  ranged  from  4.  8  x  10  to 
3.  6  x  10'2. 

These  experiments  have  demonstrated  the  application  of  the  shock 
tunnel  to  rarefied- gas  dynamic  research.  It  should  be  noted  that  the  11-  by 
15-inch  shock  tunnel,  used  in  this  program,  was  never  designed  for  low- 
density  studies.  In  fact,  it  has  a  two- stage  nozzle  with  a  rectangular 
cross-section,  which  is  relatively  poor  for  low-density  flows  with  their 
associated  thick  nozzle-wall  boundary  layers.  Nevertheless,  no  significant 
changes  in  operating  techniques  or  instrumentation  were  required  for  the 
present  program.  The  shock  tunnel  instrumentation  developed  in  recent 
years  proved  entirely  satisfactory  for  measuring  pressures  and  heat-transfer 
rates.  Based  on  the  success  of  these  experiments,  a  detailed  study  of  the 
requirements  for  a  hypersonic,  low-density  wind  tunnel  has  been  undertaken 
to  explore  additional  applications  of  the  shock  tunnel  to  rarefied  gasdynamic 
research. 

In  the  present  report,  the  various  flow  regimes  between  continuum  and 

12 

free-molecule  flow,  as  classified  by  Hayes  and  Probstein  ,  are  described. 

Some  of  the  pertinent  theoretical  solutions  found  for  these  regimes  are 

indicated.  In  addition,  a  new  theoretical  treatment  of  viscous  shock-layer 

13 

flow  by  H.  K.  Cheng  ,  which  is  used  to  correlate  the  hemisphere-cylinder 
data,  is  briefly  discussed.  Then  the  present  experimental  program  is 
described.  The  results  are  compared  with  theoretical  predictions  and,  for 
the  hemisphere  cylinder,  with  data  obtained  by  Ferri,  Zakkay,  and  Ting*® 
in  a  conventional,  blow-down  wind  tunnel.  Next,  the  results  of  the  study  of 
low-density  wind  tunnel  requirements  are  presented  indicating  that  the  shock 
tunnel  possesses  many  advantages,  some  of  them  unique,  that  make  it  well 
suited  for  research  in  hypersonic  rarefied  gasdynamics.  Finally,  future 
shock  tunnel  developments  which  will  enhance  its  capabilities  for  low-density 
research  are  described. 
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LOW -DENSITY  FLOW  PHENOMENA 


It  has  been  only  during  the  past  several  years  that  studies  of  rarefied- 

gas  flows  have  been  extended  to  hypersonic  speeds.  For  many  years  fol- 

14 

lowing  the  original  work  of  Tsien  ,  the  regimes  of  rarefied  gasdynamics 

were  classified  as  ’’continuum  flow,  M  ;,slip  flow,  ”  ’’transition  flow,  ”  and 

12 

’’free  molecule  flow.  M  Recently,  Hayes  and  Probstein  have  suggested  a 
classification  of  hypersonic  rarefied- gas  flows  as  ’’boundary-layer  regime,  " 
’’vorticity-interaction  regime,”  ’’viscous -layer  regime,”  ’’incipient  merged- 
layer  regime,”  ’’fully-merged  layer,  ”  ’’transitional  layer  regime,”  "first- 
order  collision  theory  regime,  ”  and  ”f ree-molecule  flow  regime.”  These 
various  regimes  have  been  treated  either  in  general  or  in  detail  in  References 
6,  9,  12,  15-18;  however,  the  classification  is  still  somewhat  tentative. 

As  indicated  by  the  references  cited  above,  there  is  a  growing  body  of 
literature  devoted  to  the  theory  of  hypersonic  rarefied-gas  flows.  Experimental 
results  for  low-density  flows  at  Mach  numbers  greater  than  6  are  relatively* 
rare,  however.  References  19  and  20  report  heat-transfer  measurements  to 
transverse  cylinders  at  Mach  numbers  of  about  6  and  stagnation  temperatures 
of  about  300° K.  Stagnation-point  heat-transfer  measurements  at  low  Reynolds 
numbers  on  hemisphere  cylinders  at  a  Mach  number  of  about  6  and  stagnation 
temperatures  from  4500°K  to  6000*K  have  been  reported  in  Reference  21. 
Similar  measurements  at  Mach  8  and  stagnation  temperatures  of  890*K  and 
1280°K  are  reported  in  Reference  10.  Measurements  of  heat  transfer  to  fine 
wires  in  free-molecule  flow  at  Mach  numbers  from  10  to  18  are  reported  in 
Reference  22.  These  experiments  were  made  in  a  hypersonic  shock  tunnel 
at  stagnation  temperatures  from  1800#K  to  2850*K.  The  latter  three  refer¬ 
ences  and  the  present  work  constitute  the  only  high-temperature,  hypersonic, 

low-density  data  known  to  the  authors.  The  experiments  reported  in  Refer- 
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ences  10  and  21  are  devoted  to  studying  the  vorticity  interaction  effects  , 
while  the  experiments  of  Reference  22  deal  with  free-molecule  flow  about  a 
transverse  cylinder.  The  present  study  extends  from  the  region  of  boundary- 
layer  flow  into  the  fully-merged  layer  regime. 
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Before  describing  the  present  experiments,  the  low-density  phenomena 

encountered  in  the  stagnation  region  of  a  blunt  body  at  hypersonic  speeds  will 

be  described.  This  description  will  also  very  briefly  summarize  the  analyses 

contained  in  References  6-10,  12  and  15-18*  The  regimes,  as  previously  cited, 

extend  from  continuum  flow  to  free  molecule  flow  as  the  density  is  decreased. 

13 

The  theoretical  analysis  of  Cheng  ,  which  is  presently  available  only  as  a 
C.A.L.  memorandum,  will  be  described  separately. 

Boundary-Layer  Regime 

At  high  Reynolds  numbers  viscous  effects  are  restricted  to  a  thin 

boundary  layer,  and  the  rest  of  the  flow  field  between  the  body  and  the  shock 

wave  is  inviscid.  The  classical  boundary- layer  theory  (Navier-Stokes 

23 

equations)  for  stagnation-point  flows  has  been  developed  by  Sibulkin  , 

24  25  26 

Reshotko  and  Cohen  ,  Lees  ,  Fay  and  Riddell  and  others.  Sibulkin* s 

solution  is  for  laminar,  incompressible,  low-speed  flow  about  a  body  of 

revolution.  The  comparable  solution  for  two-dimensional  flow  was  obtained 
27 

by  Squire  .  Of  the  various  compressible  boundary- layer  theories,  that  of 

26 

Fay  and  Riddell  is  the  most  complete.  Their  analysis  considers  the  case 
where  the  external  inviscid  flow  may  be  dissociated  and  includes  the  effects 
of  diffusion  and  atom  recombination  in  the  boundary  layer.  This  analysis, 
which  is  applicable  to  both  transverse  cylinders  and  bodies  of  revolution,  has 
been  used  for  comparison  of  the  present  experimental  results  at  the  higher 
Reynolds  numbers. 

Vorticity-Interaction  Regime 

As  the  density  and  Reynolds  number  are  lowered,  the  boundary- layer 
thickness  increases  until  the  vorticity  in  the  inviscid  flow  begins  to  affect 
the  boundary- layer  structure.  This  vorticity  is  generated  by  the  curved  bow 

shock  wave.  The  vorticity-interaction  effect  has  been  treated  by  Hayes  and 

12  821  10 
Probstein  ,  Hoshizaki  '  ,  and  Ferri,  Zakkay  and  Ting  .  The  first  two 

analyses  assume  incompressible  flow  between  the  shock  and  the  body,  while 

the  third  analysis  includes  compressibility  effects.  A  solution  for  this  regime 

13 

also  results  from  Cheng*s  theory  which  includes  compressibility  effects. 

In  this  regime  the  shock  wave  is  assumed  to  be  thin  with  negligible  viscous 
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and  heat  conduction  effects.  It  should  be  noted  that  for  the  transverse  cyl¬ 
inder,  the  vorticity  interaction  is  a  second-order  effect  (in  Reynolds  number) 
while  it  is  a  first-order  effect  for  the  sphere  or  hemisphere  cylinder^. 

Viscous-Layer  Regime 

At  Reynolds  numbers  less  than  those  required  for  the  validity  of  the 
vorticity-interaction  approach,  the  shock  layer  is  a  fully  viscous  contimuum 
region  governed  by  the  Navier-Stokes  equations,  and  the  shock  wave  is  a 
discontinuity  across  which  the  Rankine-Hugoniot  relations  apply.  The  latter 
condition  establishes  a  limit  to  this  regime  by  requiring  that  the  viscous 
shear  be  small  compared  to  the  tangential  momentum  transport  across  the 
shock  wave. 

Constant-density  solutions  have  been  obtained  by  Probstein  and  Kemp^, 

17  28  12 

Oguchi  *  and  Hayes  and  Probstein  .  Compressible  flow  solutions  have 
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been  reported  by  Herring  and  Ho  and  Probstein  .  These  analyses  solve  the 
Navier-Stokes  equations  using  conditions  behind  the  shock  wave  as  the 
boundary  conditions  at  the  outer  edge  of  a  boundary-layer-type  flow.  These 
analyses  are  also  restricted  to  the  immediate  neighborhood  of  the  stagnation 
point. 

Incipient  Merged  Layer  Regime 

Here  the  shock  wave  can  no  longer  be  considered  as  a  discontinuity 

obeying  the  Rankine-Hugoniot  relation,  and  the  shock-layer  flow  is  still  con- 
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sidered  as  a  continuum.  Two  possible  approaches  are  to  consider  the  shock 
as  a  discontinuity  obeying  the  conservation  laws  in  which  the  viscous  stresses 
and  heat  conduction  are  accounted  for,  or  to  use  the  complete  Navier-Stokes 
equations  to  obtain  a  solution  which  includes  the  shock-wave  structure  and 
has  free-stream  conditions  as  outer  boundary  conditions.  The  latter  approach 
is  detailed  by  Probstein  and  Kemp^,  who  present  a  numerical  solution  for 
the  sphere. 

Fully  Merged  Layer  and  Transitional  Layer  Regimes 

This  regime  is  bounded  by  the  limits  for  the  incipient  merged  layer 
regime  and  the  first-order  collision  theory  regime.  The  shock  wave  is  no 
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longer  a  discontinuity  and  the  shock  layer  is  in  either  an  almost  continuum 

regime  or  an  almost  first-order  collision  regime.  Accurate  treatment  of 

this  flow  regime  requires  the  full  formulation  of  kinetic  theory.  Failing 
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this,  Hayes  and  Probstein  suggest  that  the  best  approach,  for  engineering 
purposes,  is  to  apply  the  complete  Navier-Stokes  equations.  The  limits  of 
validity  of  this  approach  are  unknown. 

First-Order  Collision  Theory  Regime 

In  this  regime  the  Knudsen  number,  appropriately  defined,  is  large, 

but  not  large  enough  for  free-molecule  flow  concepts  to  be  valid.  This  flow 

region  can  be  analyzed  by  considering  only  collisions  between  free-stream 

molecules  and  molecules  re- emitted  from  the  surface.  Such  analyses  have 
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been  made  by  Lunc  and  Lubonski  ,  Baker  and  Charwat  ,  Hammerling  and 

31*32  33 

Kivel  ,  Liu  ,  and  Willis  . 

Free-Molecule  Flow  Regime 

Here  all  interm olecular  collisions  are  neglected  and  only  the  interaction 

of  free-stream  molecules  with  the  surface  is  considered.  There  is  a  large 

body  of  literature  on  this  subject.  Recent  general  treatments  of  this  flow 

,34 

regime  are  given,  for  example,  by  Schaaf  and  Chambre  ,  Schaaf  and 
35  i  7 

Talbot  ,  and  Hayes  and  Probstein 

Low- Density  Regime  Boundaries 

The  limits  of  validity  of  these  various  regimes  have  been  discussed  in 

References  6,  12,  15,  16  and  18.  For  flow  in  the  stagnation-point  region  of 

18 

a  blunt  body  they  are  given  by  Probstein  as  follows: 


Boundary  Layer 

Vorticity  Interaction 

Xs/R6  <<  £ 2 

Viscous  Layer 

xs/*b«e3/ 

Incipient  Merged  Layer 

Fully  Merged  Layer 

*  s/*k  <  £ 

Transitional  Layer 

^s/ ^b  ^  £ 
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First-Order  Collision  Theory 


a  'I'm?  ’  rr  7 

Free-Molecule  Theory  >>  — /-r— 

▼  1  *w 

where  ^  is  the  body  radius,  £  is  the  density  ratio  across  the  shock  wave, 
Poo/fi$  *  and  A*©  ,  A$  and  A^  are  the  mean  free  paths  in  the  free-stream 
flow,  behind  the  shock  wave,  and  at  the  body,  respectively.  These  regimes 
are  shown  in  Figure  1  which  is  taken  from  Reference  18. 

On  the  basis  of  these  boundaries,  the  present  experiments  extend 
slightly  into  the  fully  merged  layer  regime  for  the  transverse  cylinder  and 
up  to  the  boundary  between  the  incipient  merged  and  fully  merged  layers  for 
the  hemisphere  cylinder. 

* 

Low-Density  Shock-Layer  Theory  of  Cheng 
13 

H.  K.  Cheng  has  recently  developed  an  approach  to  the  hypersonic 
blunt-body  problem  utilizing  the  shock- layer  concept  and  including  viscous 
and  heat- conduction  effects  in  the  shock  layer  as  well  as  in  the  conservation 
relations  of  the  shock  itself.  The  work  is  as  yet  published  only  in  Reference 
13  which  has  a  limited  distribution.  This  investigation  deals  with  flows  in 
which  the  density  is  so  low  that  the  transport  process  may  be  as  important  as 
the  convective  process,  but  not  so  low  that  the  continuum  concept  becomes 
invalid.  Such  an  approach  should  be  valid  until  the  smallness  of  the  density 
requires  a  kinetic  theory  treatment  of  the  gas.  The  validity  of  the  approach 
is  further  strengthened  by  the  fact  that  the  theory  aggrees  with  free-molecule 
flow  theory  in  the  low-density  limit.  At  present  only  the  stagnation  region 
for  a  body  of  revolution  has  been  analyzed.  It  is  planned  to  extend  the  study 
to  two-dimensional  flows. 

Using  this  approach,  Cheng  has  found  that  the  problem  may  be  sepa¬ 
rated  into  two  regimes.  In  one  regime,  where  the  Rankine-Hugoniot  shock 

12 

relations  hold,  the  vorticity-interaction  theory  of  Hayes  and  Probstein 

♦  "  — 

This  work  has  been  performed  under  the  sponsorship  of  the  U.  S.  Navy 
Office  of  Naval  Research  (Contract  Nonr-2653(0Q) )  and  the  U.  S.  Air  Force 
Rome  Air  Development  Center  (Contract  AF  30(602)-2267).  The  presentation 
here  has  been  taken  from  Reference  13. 
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based  on  the  boundary-layer  approximation  is  found  to  be  adequate  even  in 
the  range  of  validity  of  viscous-layer  theory.  For  the  other  regime,  where 
the  viscous  and  heat-conduction  effects  on  the  shock  relations  cannot  be 
neglected,  an  analytic  solution  of  the  viscous  shock  layer  has  been  obtained. 
This  solution  is  subject  to  assumptions  regarding  the  viscosity-temperature 
variation,  the  magnitude  of  €*=  fis / 2/)$  Hs  ,  and  the  order  of  K* * 
where  C  is  the  Chapman-Rubesin  viscosity-temperature  coefficient^.  Also, 
a  perfect  gas  with  constant  specific  heats  has  been  assumed.  This  solution 
is  valid  in  the  fully-merged  layer  regime. 

37  38 

As  with  the  inviscid  shock  layer  theory  *  ,  Cheng*s  analysis  of  the 

viscous  shock  layer  requires  a  high  shock  compression  ratio,  which  in  turn 
requires  a  strong  shock  and  a  specific  heat  ratio  not  far  from  unity,  i.  e.  , 

€ '  <<  1  .  The  flow  model,  because  of  the  simplicity  resulting  from  the 

shock-layer  concept,  permits  reduction  of  the  problem  to  one  comparable 
to  that  of  the  boundary- layer  approximation.  However,  the  formulation, 
which  has  been  developed  to  within  an  error  of  order  € '  ,  differs  from  the 
boundary-layer  theory  in  that  the  normal  pressure  gradient  is  not  generally 
negligible  and  in  that  the  outer  edge  of  the  shock  layer  lies  at  a  finite  distance 
from  the  body  surface.  Implicit  in  the  flow  model  is  the  continuum  assumption 
since  the  compressible,  Navier-Stokes  equations  are  assumed  to  govern  the 
flow. 

The  solutions  obtained  by  Cheng  will  be  compared  with  the  present  hemi¬ 
sphere-cylinder  data  as  well  as  the  data  reported  by  Ferri  et  al^. 
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EXPERIMENTS 


Hypersonic  Shock  Tunnel 

4 

The  C.A.L.  11-  by  15-inch  hypersonic  shock  tunnel  consists  of  a 

high-pressure  driver  tube  3.5  inches  I.  D.  by  14  feet  long,  a  driven  tube 

3  inches  I.  D.  by  28  feet  long,  and  a  three- stage  hypersonic  nozzle  (Figure  2) 

having  a  variable  entrance  area  in  the  final  stage.  The  driver  and  driven 

tubes  are  designed  for  30,  000  psi  operation  using  either  a  single  or  double 

diaphragm  between  the  two  sections.  The  driven  tube  has  instrumentation 

ports  located  at  intervals  along  it  as  well  as  charging  and  charge  pressure 

measuring  ports.  The  driven  tube  is  separated  from  the  nozzle  system  by 

a  Mylar  diaphragm.  Diaphragm  thicknesses  from  1/2  to  2  mils  are  used 

depending  upon  the  initial  driven  tube  charge  pressure.  The  purpose  of 

this  diaphragm  is  to  separate  the  two  sections  in  order  to  maintain  the 

proper  driven  tube  pressure  while  permitting  evacuation  of  the  test  section 

39 

to  promote  efficient  tunnel  starting 

The  first  nozzle  stage  consists  of  a  convergent -diver gent  nozzle  with 
a  contraction  ratio  of  19:1-  This  is  designed  so  as  to  completely  reflect  the 
incident  shock  wave,  thus  producing  a  region  of  nearly  stagnant  air  at  high 
temperature  and  pressure  behind  the  reflected  shock.  The  divergent  portion 
is  contoured  to  produce  uniform  flow  at  about  Mach  4.  The  second  stage 
consists  of  a  wedge-shaped  plate  placed  at  an  angle  of  attack  of  10*  to  the 
flow  emerging  from  the  first-stage  nozzle.  A  Prandtl-Meyer  expansion  from 
the  wedge  deflects  the  air  flow.  Since  the  diaphragm  particles,  due  to  their 
higher  momentum,  cannot  turn  with  the  air  stream,  they  are  separated  from 
the  flow.  The  third  stage  is  placed  at  10*  from  the  axis  of  the  first  stage  to 
accept  the  centrifuged  flow.  This  stage  consists  of  a  straight-walled,  15- 
degree  included-angle  nozzle  terminating  in  an  11-  by  15-inch  test  section. 
The  test-section  Mach  number  can  be  varied  by  changing  the  entrance  area 
of  the  third- stage  nozzle.  All  tests  in  the  present  investigation  were  con¬ 
ducted  using  one  fixed  area  ratio. 
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Operational  Procedure 

Prior  to  shock  tunnel  operation,  the  test  section  is  evacuated  to 

approximately  1  micron  of  mercury*  At  the  same  time  the  driven  and 

driver  tubes  are  evacuated  to  prevent  possible  contamination  of  the  charge 

gases.  The  driven  tube  is  then  pressurized  with  the  gas  chosen  for  the 

test  (in  these  tests  dry-air  was  used).  Finally,  the  driver  is  pressurized 

with  either  hydrogen  or  helium  to  a  pressure  which  permits  the  proper 

4 

driver  to  driven  pressure  ratio  for  "tailored  interface"  operation  . 

Two  methods  of  test  initiation  were  used  during  this  series.  For  the 
transverse-cylinder  experiments,  a  single  diaphragm  was  used  between 
driver  and  driven  tubes.  After  the  driven  pressure  was  set,  the  driver  was 
pressurized  until  the  diaphragm  burst.  -The  burst  pressure  of  the  metal 
diaphragm  was  controlled  by  the  commonly-used  technique  of  pressing  an 
X  in  the  diaphragm  to  a  specific  depth,  At  low  pressures,  when  Mylar  was 
used  as  the  main  diaphragm,  the  material  strength  was  the  determining 
factor  for  burst  pressure.  This  practice  gave  burst  pressure  ratios  which 
varied  from  run  to  run. 

In  order  to  obtain  predictable  and  consistent  burst  pressures  for  the 

* 

later  hemisphere-cylinder  tests,  the  double  diaphragm  technique  was  tried. 
Two  diaphragms  of  equal  strength  were  placed  between  the  driver  and  driven 
tubes  such  that  a  small  volume  of  gas  at  about  half  the  driver  pressure  could 
be  held  between  them.  The  diaphragms  were  scribed  so  that  they  would,  if 
used  singly,  burst  at  approximately  60  to  70  percent  of  the  expected  driver 
pressure.  By  pressurizing  the  region  between  the  diaphragms  to  one-half  the 
driver  pressure  and  then  loading  the  driver,  the  exact  pressure  ratio  desired 
between  driver  and  driven  tubes  could  be  obtained.  To  initiate  the  run,  the 
pressure  between  the  two  diaphragms  was  reduced.  This  procedure  was  used 
successfully  for  driver  pressures  as  low  as  50  psig.  Diaphragm  burst  pres¬ 
sures  were  extremely  repeatable  and  could  easily  be  varied. 

n - 

The  first  successful  application  of  this  technique  brought  to  the  authors' 
attention  was  by  R.  N,  Cox  of  the  A.  R.  D.  E.  at  Fort  Halstead,  England.  The 
first  application  at  C.A.L.  was  in  the  24"  hypersonic  shock  tunnel  operated 
by  the  Experimental  Facilities  Division. 
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Rupture  of  the  diaphragm  between  the  driver  and  driven  tubes  results 
in  the  formation  of  a  shock  wave  which  travels  through  the  driven  gas 
heating  and  compressing  it.  The  shock  wave  reflects  from  the  end  of  the 
driven  tube  and  again  passes  through  the  test  air  further  increasing  the 
temperature  and  pressure.  The  resulting  pressure  causes  the  Mylar 
diaphragm  between  the  driven  tube,  and  nozzle  to  rupture,  permitting  the 
heated  test  air  to  pass  through  the  three-stage  hypersonic  nozzle. 

4 

By  utilizing  the  tailored-interface  technique  developed  at  C.  A.  L.  , 
steady  flow  may  be  maintained  in  the  test  section  for  several  milliseconds. 

With  room-temperature  helium  as  the  driver  gas,  tailored-interface  opera¬ 
tion  will  produce  stagnation  temperatures  in  the  range  of  1800°K  to  2000°K 
for  a  test  duration  of  about  6-1/2  milliseconds.  Using  room-temperature 
hydrogen  as  the  driver  gas,  stagnation  temperatures  from  3600°K  to  4000°K 
are  produced  for  flow  periods  of  about  3-1/2  milliseconds. 

Flow  Conditions 

The  parameters  recorded  for  each  test  are  the  following:  initial  gas 
pressure  in  the  driven  tube;  driver  gas  pressure  at  the  time  the  test  is 
initiated;  pressure  behind  the  incident  and  reflected  shock  wave;  incident 
shock  wave  speed;  static  pressure  on  the  test  section  sidewall;  the  stagnation- 
point  heat-transfer  rate;  the  stagnation- point  temperature  change;  test  section 
pressure  reading  before  vacuum  system  is  closed;  leak  rate  of  test  section; 
time  from  closing  vacuum  system  to  test  initiation. 

The  pressure,  temperature,  enthalpy,  and  entropy  of  the  air  behind 
the  reflected  shock  wave  (which  are  the  stagnation  or  reservoir  conditions 
for  the  nozzle)  are  calculated  from  the  measured  shock  speed  and  the  initial 
air  pressure  using  tables  and  charts  for  real  air  in  thermodynamic  equilib¬ 
rium,  for  example,  References  40-42.  The  agreement  of  the  calculated  and 
measured  stagnation  pressure  provides  a  check  as  to  the  validity  of  the  other 
calculated  parameters.  The  test  section  Mach  number  is  calculated  using 
the  measured  free-stream  static  pressure  and  the  stagnation  pressure 
assuming  isentropic  flow  in  the  nozzle.  In  most  tests,  the  measured  stagna¬ 
tion  pressure  and  the  value  of  stagnation  pressure  calculated  using  the  incident 
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shock  Mach  number  are  not  in  agreement.  Although  the  disagreement  is 
appreciable,  Figure  3,  the  effect  on  the  reduced  test-section  conditions  and 
data  is  not  great.  Figure  4  is  a  plot  of  the  free- stream  Mach  number  calcu¬ 
lated  using  the  static  pressure  and  the  stagnation  or  reservoir  pressure 
both  as  measured  and  as  calculated.  Several  possible  reasons  for  the 
pressure  difference  exist.  Inaccurate  calibration  of  piezo-electric  crystal 
transducers  whose  output  signal  is  low  at  these  pressure  ranges  does  con¬ 
tribute  some  error.  However,  it  is  felt  that  this  does  not  account  for  all  of 
the  error.  The  loss  of  mass  flow  through  the  convergent -diver gent  nozzle 
was  studied  to  approximate  that  contribution.  In  the  flow  time  in  question 
there  appeared  to  be  little  or  no  additional  error.  A  small  contribution  is 
possible  if  shock  attenuation  is  not  accounted  for  between  the  final  wave 
speed  transducer  and  the  end  of  the  shock  tube.  This  distance,  however, 
is  so  small  (5  inches)  that  even  if  attentuation  were  not  considered  here, 
the  effect  upon  the  value  of  the  incident  shock  Mach  number  used  in  finding 
the  stagnation  pressure  ratio  would  be  negligible.  One  further  possibility 
considered  is  that  waves  are  initiated  by  diaphragm  flaps  which  rebound 
from  the  driven  tube  walls  after  the  diaphragm  ruptures.  This  effect  could 
vary  depending  upon  the  mass  of  the  diaphragm  flap.  It  would  appear  that 
any  waves  produced  would  have  to  lag  behind  the  incident  shock  since  no 
noticeable  effect  on  wave  speed  has  been  witnessed. 

Figure  3  indicates  that  a  consistent  difference  in  the  measured  and 
calculated  pressures  occurred  both  in  the  single-diaphragm  and  double¬ 
diaphragm  tests.  It  is  significant  that  the  measured  pressure  was  always 
less  than  the  calculated  pressure  for  the  single- diaphragm  runs  and  always 
greater  for  the  double- diaphragm  tests.  This  seems  to  indicate  clearly 
that  conditions  at  the  diaphragm  station  are  related  to  the  discrepancy. 
Whether  the  difference  is  related  to  diaphragm  opening  time  or  by  constric¬ 
tion  of  the  flow  by  rebounding  diaphragm  flaps  has  not  been  determined. 

This  matter  requires  further  investigation.  In  any  case,  a  consideration 
of  the  possible  causes  of  the  pressure  difference  seems  to  indicate  that  the 
measured  pressure  be  used  as  the  flow  stagnation  pressure. 
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This  discrepancy  between  the  measured  and  calculated  pressures  was 
only  one  of  the  difficulties  encountered  in  making  low  pressure  shock  tunnel 
runs.  In  addition,  there  were  such  matters  as  increased  shock  wave  attenu¬ 
ation  (although  it  was  still  small  enough  to  be  unimportant),  the  difficulty  in 
obtaining  satisfactory  low-pressure  diaphragms  (this  was  largely  overcome 
by  use  of  the  double- diaphragm  technique),  and  the  thick  nozzle- wall  boundary 
layers  (particularly  in  the  first-stage  nozzle  which  has  a  1.  5-  by^  12.  25-inch 
exit  cross-section).  In  spite  of  these  problems,  hypersonic  low-density 
flow  was  obtained  which  was  satisfactory  for  stagnation-point  heat-transfer 
experiments. 

As  an  indication  of  the  degree  of  rarefaction  which  has  been  attained 
during  this  series  of  runs  a  curve  of  Knudsen  number  versus  Reynolds 
number  behind  the  shock  wave,  Res  ,  has  been  plotted  in  Figure  5.  In 
calculating  the  Knudsen  number,  the  mean  freepath,  Xs  ,  behind  the  normal 
shock  and  body  radius,  R^  ,  have  been  used.  The  various  rarefied  flow 
regimes  are  indicated  on  this  figure. 

For  the  experiments  in  which  helium  was  used  as  the  driver  gas,  the 
range  of  flow  conditions  is  given  in  Table  I: 


TABLE  I 

FLOW  CONDITIONS  FOR  HELIUM  DRIVER  GAS  EXPERIMENTS 


RESERVOIR  PRESSURE,  f0 
RESERVOIR  TEMPERATURE,  T0 
FREE- STREAM  VELOCITY, 

FREE-STREAM  MACH  NUMBER, 

FREE-STREAM  UNIT  REYNOLDS  HUMBER,  Re^/L 
REYNOLDS  NUMBER  BEHIND  SHOCK ,  Res 
KNUDSEN  NUMBER,  Xs/Rk 


2.4-58.0  ATM. 
1500-2200  CK 
5500-6700  FT/SEC 
9.7-11.6 

1460-21,300  PER  INCH 
1 1 - 1 000 
0.03 6-. 00048 


Foi  the  tests  using  hydrogen  as  the  driver  gas,  the  range  of  flow 
conditions  is  listed  in  Table  II: 
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TABLF.  II 

FLOW  CONDITIONS  FOR  HYDROGEN  DRIVER  GAS  EXPERIMENTS 


RESERVOIR  PRESSURE,  p0 
RESERVOIR  TEMPERATURE,  T0 
FREE-STREAM  VELOCITY;  U 
FREE-STREAM  MACH  NUMBER,  M«, 
FREE-STREAM  UNIT  REYNOLDS  NUMBER, 
REYNOLDS  NUMBER  BEHIND  SHOCK,  R»$ 
KNUDSEN  NUMBER,  A S/R b 


3.48-58.1  ATM. 
3200-3700  °K 
9000-10,000  FT/ SEC 
8.4-10.3 

561-3410  PER  INCH 
16-86 

0.0I2-.002 


Free- stream  Mach  number  was  calculated  by  isentropically  expanding 
the  stagnation  conditions  to  the  measured  free- stream  static  pressure.  A 
correction  of  the  static  pressure  was  made  to  account  for  the  test  section 
leak  rate  after  the  vacuum  system  was  closed.  The  leak  rate  was  noted  and 
the  time  between  closing  of  the  vacuum  system  and  rupture  of  the  diaphragm 
was  used  to  compute  the  back  pressure  on  the  static  pressure  transducer. 

A  check  upon  this  method  of  obtaining  the  free- stream  Mach  number  was 
made  using  the  results  from  a  pitot  probe  placed  in  the  air  flow  of  the  test 
region.  The  Mach  number  based  on  the  ratio  of  pitot-to-reservoir  pressures 
agreed  with  that  calculated  from  the  static  pressure  within  the  scatter  of  the 
results,  Figure  4. 

A  certain  amount  of  variation  in  reservoir  temperatures  can  be  seen 
in  the  tables  (I  and  II).  Obviously  this  variation  also  appears  in  the  stagna¬ 
tion  heat  transfer,  f  .  It  results  from  the  initial  inconsistency  of  the  pres¬ 
sure  ratio  at  which  the  diaphragm  between  driver  and  driven  tubes  burst. 
Thus,  a  variation  in  shock  Mach  number  occurs  from  run  to  run  which 
produces  a  variation  in  temperature  behind  the  reflected  shock  wave. 
Primarily,  this  inconsistency  occurred  during  the  transverse  cylinder 
tests.  The  later  hemisphere-cylinder  tests  utilized  the  double  diaphragm 
technique  which  produced  more  repeatable  flow  conditions. 

Tables  I  and  II  also  indicate  a  variation  of  Reynolds  number  and  free- 
stream  Mach  number.  The  variation  of  Mach  number  with  Reynolds  number 
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is  caused  by  the  influence  of  the  nozzle-wall  boundary  layer.  This  is 
illustrated  in  Figure  4  which  shows  test  section  Mach  number  as  a  function 
of  nozzle  stagnation  pressure. 

Models  and  Instrumentation 

Two  series  of  shock  tunnel  tests  have  been  mentioned  previously. 

The  first  series  included  tests  using  three  cylindrical  models  having  nominal 

diameters  of  1/8,  1/4,  and  1  inch  (Figure  6).  The  models  were  mounted 

transversely  to  the  air  stream  with  their  axes  in  the  plane  of  symmetry  of 

'  43 

the  wedge-shaped  final  nozzle.  Thin-film  resistance  thermometers  were 
mounted  along  the  stagnation  line  of  the  cylinders.  The  1 /4-inch  and  1-inch 
models  had  three  resistance  thermometers  mounted  about  one  inch  apart. 

The  1/8-inch  model  had  a  single  resistance  thermometer. 

For  the  second  series  a  1/2-inch  diameter  hemisphere-cylinder  with 
a  resistance  thermometer  mounted  at  the  stagnation  point  of  the  hemispherical 
portion  of  the  model  was  used.  Figure  7.  A  technique  described  in  Refer¬ 
ences  43  and  44  was  used  whereby  a  resistance  thermometer  is  formed  by 
painting  a  0.  1  micron-thick  platinum  film  upon  a  specified  backing  material, 
then  curing  the  film  by  baking.  The  thi^n-film  resistance  thermometer  was 
a  rectangle  measuring  1  / 1 6  inch  either  side  of  the  stagnation  point  by  approxi¬ 
mately  1/32  inch  wide.  A  colloidal- silver  paint  strip  which  is  baked  on  the 
model  acts  as  a  lead  to  the  wires  at  the  rear  of  the  hemisphere- cylinder. 
Resistance  thermometers,  when  used  for  quantitative  results,  must  be 
backed  by  a  dielectric  material  with  known,  constant  properties.  Conse¬ 
quently,  Pyrex  glass  was  used  in  making  all  models  tested  in  the  two  series. 
Since  the  heat  capacity  of  the  thin  platinum  film  is  negligible,  it  is  assumed 
that  the  film  temperature  is  the  instantaneous  surface  temperature  of  the 
glass;  thus  heat  lag  due  to  film  thickness  may  be  neglected. 

Calibration  of  the  thin-film  resistance  thermometer  is  done  using  a 
controlled-temperature  water  bath.  The  model,  enclosed  in  a  watertight 
bag,  is  allowed  to  soak  in  the  bath  until  it  has  reached  the  temperature  of 
the  bath.  This  is  detected  by  a  thermocouple  on  the  model  surface.  As  the 
bath  temperature  is  increased,  a  linear  rise  in  resistance  is  obtained.  The 
normal  calibration  range  is  from  room  temperature  to  50#C. 


15 


During  a  shock  tunnel  test  the  thermometers  were  operated  at  a  constant 
electrical  current  so  that  surface  temperature  changes,  which  produced 
changes  in  film  resistance,  would  result  in  a  directly  proportional  voltage 
signal.  This  signal  can  be  recorded  directly  on  an  oscillograph  from  which 
a  temperature  change  can  be  obtained  for  later  use  in  calculating  the  heat 
flux.  The  signal  can  also  be  fed  into  an  analogue  computer  network  such  as 
described  in  Reference  45.  The  network  converts  the  output  of  the  thin-film 
resistance  thermometer  to  an  electrical  signal  proportional  to  the  surface 
heat  flux.  An  oscillograph  record  of  heat-flux  data  is  available  immediately 
after  a  test  run  using  the  analogue  network. 
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EXPERIMENTAL  RESULTS 


Transverse-Cylinder  Tests 

In  the  data  reduction  and  presentation,  the  measured  stagnation-point 

heat-transfer  rate,  £  ,  has  been  nondimensionalized  with  respect  to  the 

quantity  This  is  equivalent  to  the  conventional  Stanton  number 

for  the  case  of  a  cold  wall,  <<  T0  .  The  data  are  presented  in  Figure  8 

as  the  nondimensionalized  heat  transfer  multiplied  by  the  Reynolds  number 

based  on  flow  conditions  behind  the  shock  wave  and  body  radius,  j  , 

and  are  plotted  as  a  function  of  the  Reynolds  number  behind  the  shock  wave, 

* 

*  The  data  represented  by  the  open  symbols  are  experiments  for  which 

the  density  ratio  across  the  shock  wave,  €  5  fi**//*$  9  ^es  between  0.  148  and 

0.  152;  these  are  helium-driver- gas  experiments  for  which  the  stagnation 

temperature  is  in  the  range  from  1800°K  to  2100°K.  At  the  higher  Reynolds 

numbers  the  data  are  seen  to  be  in  fairly  good  agreement  with  the  boundary- 

26  ^ 

layer  theory  of  Fay  and  Riddell  as  calculated  for  these  flow  conditions. 

At  Reynolds  numbers  below  30  the  data  tend  to  fall  increasingly  below  the 
boundary- layer  theoretical  prediction. 


In  the  very  low  Reynolds  number  region  of  these  experiments  it  was 
difficult  to  obtain  precisely  repeatable  test  conditions,  because  of  the  diffi¬ 
culty  involved  in  making  thin  diaphragms  that  would  burst  at  precisely  the 
desired  pressure.  Consequently,  some  of  the  data  fall  outside  the  density 
ratio  range  from  0.  148  to  0.  152,  but  rather  were  in  the  range  from  0.  157  to 
0.  163.  These  data  are  indicated  by  the  half-filled  symbols. 


The  hydrogen-driver- gas  experiments  are  shown  as  the  solid  symbols, 
Figure  8.  Only  the  1 /4-inch  model  has  been  tested  at  these  conditions.  For 


Utilizing  the  continuity  equation  across  a  normal  shock  wave,  fl+Um  mPsUgf 
the  Reynolds  number  behind  the  shock  wave  can  be  written  as  ^ 

Subsequent  to  these  experiments,  the  double- diaphragm  operating  technique 
was  introduced  for  the  hemisphere- cylinder  tests. 
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these  experiments  the  density  ratio  was  between  0.  10  and  0.  1Z  for  stagnation 

temperatures  from  3200°K  to  3700°K.  The  boundary- layer  theory  has  been 

calculated  for  a  density  ratio  of  0.  11.  The  data  are  seen  to  lie  somewhat 

above  the  theory  for  this  series  of  experiments,  whereas  for  the  lower 

stagnation  temperature  tests,  the  data  fell  below  the  theory  in  the  same 

Reynolds  number  range.  In  view  of  the  low  stagnation  pressure  levels  of 

these  experiments,  it  was  considered  possible  that  the  air  flow  might  not  be 

in  thermodynamic  equilibrium.  The  effects  of  freezing  of  the  flow  on  the 

test- section  conditions  and  on  the  stagnation-point  heat-transfer  rate  have 

46 

been  investigated  by  Bray  .  He  has  shown  that  freezing  in  the  nozzle 
causes  a  reduction  in  the  measured  heat-transfer  rate.  At  temperatures  of 
3700°K  and  less,  nonequilibrium  effects  would  be  expected  to  be  small.  The 
agreement  of  the  present  data  with  theory  indeed  indicates  that  these  effects 
were  apparently  negligible,  and  the  assumption  that  the  flow  is  in  equilibrium 
is  reasonably  accurate. 

It  is  observed  in  Figure  8  that  the  low  stagnation  temperature  data 
(open  and  half-filled  symbols)  tend  to  lie  below  the  boundary-layer  theory 
for  all  Reynolds  numbers  below  30,  whereas  the  high-temperature  data 
(solid  symbols)  are  in  good  agreement  with  theory  at  Reynolds  numbers 
down  to  16.  This  results  from  the  difference  in  the  degree  of  rarefaction 
of  the  two  flows  at  a  given  Reynolds  number.  For  example,  at  a  Reynolds 
number  of  20,  the  low-temperature  flow  has  a  Knudsen  number  of  0.  02, 
while  the  high-tempe rature  flow  has  a  Knudsen  number  of  only  0.  01.  Hence, 
the  one  flow  is  almost  twice  as  rarefied  as  the  other  at  a  given  Reynolds 
number. 

For  the  sake  of  comparison  at  the  lowest  Reynolds  numbers,  heat- 

transfer  rates  have  been  calculated  according  to  classical  free-molecule 
12 

theory  .  These  are  shown  for  accommodation  coefficients  of  1,  1/2,  and 
1/4  (Figure  8).  Although  the  lower  stagnation  temperature  data  (open  and 
half-filled  symbols)  appear  to  tend  toward  agreement  with  free-molecule 
theory  for  an  accommodation  coefficient  of  about  1/4,  it  would  be  erroneous 
to  conclude  that  the  data  indicate  an  accommodation  coefficient  substantially 
below  unity.  Even  for  the  lowest  Reynolds  number  experiments,  the  Knudsen 
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number,  which  is  about  0.  02,  is  about  two  orders  of  magnitude  below  that  of 
f ree-molecule  flow,  Kn  >1.  At  the  lowest  Reynolds  numbers  (Figure  5),  the 
flow  is  in  the  fully-merged  region  as  described  in  References  6,  12  and  18. 

It  should  be  noted  that  in  the  Region  of  Knudsen  numbers  above  0.  01,  the  use 
of  the  inviscid  Rankine-Hugoniot  equations  to  calculate  flow  conditions  behind 
the  shock  wave  is  highly  questionable. 

The  various  symbols  (O^D)  used  in  Figure  8  indicate  the  three  heat- 
transfer  gauges  spaced  about  1  inch  apart  along  the  span  of  the  1/4-  and  1-inch 
diameter  models,  Figure  6.  The  purpose  of  the  three  gauges  was  to  check  the 
spanwise  flow  uniformity.  No  systematic  variation  greater  than  the  data 
scatter  can  be  detected  in  Figure  8. 

Presentation  of  stagnation-point  heat  transfer  in  the  form  of  Figure  8 

shows  a  great  sensitivity  to  the  density  ratio,  £  =  flu/Pt .  The  recent  work 
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by  Cheng  has  shown  that  much  of  this  sensitivity  can  be  removed  by  pre¬ 
senting  the  information  in  the  form  of  Figure  9.  Here  the  stagnation-point 
heat  transfer  is  presented  as  the  heat-transfer  coefficient 


°H  =  U) 

plotted  versus  the  product  of  the  density  ratio  and  the  Reynolds  number  behind 
the  shock,  €R*S  .  Note  that  the  high  stagnation  temperature  data  (solid  sym¬ 
bols)  and  the  low- temperature  data  (open  symbols)  and  the  corresponding 
boundary- layer  theories  are  in  close  agreement  when  presented  in  this  manner. 
As  before,  this  presentation  does  not  reveal  the  difference  in  degree  of  rare¬ 
faction  of  the  flow  for  these  two  series  of  experiments.  Hence,  the  low- 
temperature  data  tend  to  depart  from  the  boundary- layer  theory  below  values 
of  €  R*s  of  4,  whereas  the  high- temperature  data  are  still  in  good  agree¬ 
ment  with  this  theory  at  €  Rc$  =  2. 


Hemisphere- Cylinder  Tests 


The  stagnation-point  heat-transfer  data  for  the  hemisphere-cylinder 
experiments*  is  presented  in  Figure  10.  The  heat-transfer  rate  ^  is  plotted 

These  tests  were  supported  by  the  U.  S.  Navy  Office  of  Naval  Research  under 
Contract  Nonr-2653(00). 
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as  the  heat  transfer  parameter,  CH  ,  versus  the  parameter  <4  ,  which 

results  from  Cheng's  analysis  and  is  given  by 


4  = 


(2) 


where 


1^2  _  P90  %  b 

pa  Uoo  C 


(3) 


T0  is  a  reference  temperature  obtained  from  the  analysis.  For  the  present 
test  conditions,  the  corresponding  values  of  K  *  and  €'Pes  are  indicated 
along  the  abscissa  of  Figure  10.  For  strong  shock  waves  where 
and  Uc  «  ,  the  parameter  €*  «  —  -  77 —  and  the  density  ratio  across 

the  shock  wave,  €  * »  are  essentially  identical.  However,  at  con¬ 
ditions  of  the  present  experiments  these  parameters  differ  slightly.  For 
example,  in  a  typical  case,  6  =  0.  15  where  e'  *  ~2/> ~ 

The  theoretical  curves  shown  in  Figure  10  represent;  (1)  a  complete 

numerical  solution  of  the  viscous  shock- layer  equations  by  Cheng  for 

€ '  =  0.  10  and  0.  15,  (2)  Cheng's  analytic  solution  for  the  fully  viscous  flow, 

(3)  Cheng's  numerical  boundary-layer  solution  including  vorticity  (this  is 

$  • 

coincident  with  the  complete  numerical  solution  )  for  ^ €  =  0.  10  and  0.  15, 

26 

and  (4)  the  boundary- layer  theory  of  Fay  and  Riddell  which  does  not  include 
vorticity  effects.  The  analytic  solution  of  the  fully  viscous  flow  is  a  zero- 
order  result.  Higher  order  solutions  can  be  obtained  from  Cheng's  analysis 
which  will  exhibit  closer  agreement  with  the  numerical  solution  at  higher 
values  of  the  parameter,  A. 


13 


The  boundary- layer  solution  differs  from  the  complete  numerical  solution 
in  that  the  boundary- layer  approximations  are  employed  to  drop  some  terms 
in  the  Navier-Stokes  equations.  The  two  results  are  essentially  identical 
numerically  at  4  >  5. 
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The  data  are  seen  to  be  in  fairly  good  agreement  with  the  theory, 
although  the  scatter  is  greater  at  the  lower  Reynolds  number,  Figure  10. 
Although  these  experiments  extend  over  a  rather  limited  range  (only  about 
one  and  a  half  orders  of  magnitude  in  Knudsen  number  or  Reynolds  number), 
they  cover  an  important  region  from  the  vorticity-interaction  regime  through 
the  incipient  merged-layer  regime,  Figure  5.  It  is  noted  that  even  at  the 
highest  Reynolds  number  of  these  experiments,  the  vorticity-interaction 
effect  amounts  to  a  10  percent  increase  in  heat  transfer. 

Both  the  present  data  and  some  of  the  data  of  Ferri  et  al^  are  com¬ 
pared  with  Cheng's  theory  in  Figure  11.  The  theoretical  curves  have  been 
calculated  for  values  of  €*  of  0.  10  and  0.  15,  a  Prandtl  number  of  0.  71  and 
a  wall-to- stagnation  temperature  ratio,  Tw/to  ,  of  zero.  These  values  of 
€  *  bracket  the  experimental  data  which  fall  in  the  range  0.  120  £  € '  £  0.  131. 

The  wall-to- stagnation  temperature  ratios  for  the  experiments  ranged  from 

0.  150  to  0.  175  for  the  C.  A.  L.  experiments  and  had  values  of  0.  34  and  0.  46 

* 

for  the  PIBAL  experiments  .  Ferri's  data  is  seen  to  have  less  scatter  and 
to  be  in  better  agreement  with  the  theory  than  the  present  data. 

It  should  be  indicated,  that  the  PIBAL  data  is  presented  in  Reference  10 
only  as  the  ratio  of  measured  heat-transfer  rate  with  vorticity  interaction  to 
measured  heat-transfer  rate  without  voticity  interaction.  This  results  from 
the  experimental  technique  which  consists  of  testing  simultaneously  several 
models  of  different  radii  at  very  high  Reynolds  numbers,  where  the  vorticity 
interaction  should  be  zero,  and  using  the  value  of  obtained  in  these 

experiments  as  the  zero-vorticity  reference.  The  fact  that  */*r  does 
not  change  from  model  to  model  verifies  that  the  vorticity  interaction  was 
zero  for  these  experiments.  Then  the  various  sized  models  were  tested  at 
a  range  of  Reynolds  numbers, and  the  ratio  of  the  quantity  for  each 

model  to  the  value  of  for  the  largest  model  is  determined.  This  is  taken 

as  the  ratio  of  heat  transfer  with  vorticity  to  heat  transfer  without  vorticity  after 


These  values  are  inferred  from  the  fact  that  the  theory  in  Reference  10  is 
calculated  for  values  of  ° f  34  and  0.46  (see  Figure  4,  Reference 

10).  For  the  flows  conditions  of  Reference  10,  Tw/T0  can  be  taken  equal 
to  Hw/ H0  . 
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the  value  of  for  the  largest  model  is  corrected  for  vorticity  inter¬ 

action  by  ratioing  to  the  value  obtained  from  the  high  Reynolds  number  tests, 
and  the  ratio  of  the  model  radii  is  evaluated.  Thus,  Ferri  et  al^  present 
the  data  as  J  vort / Vvort  *  o  as  a  faction  of  Reynolds  number, 

Re  =■  pos  )Hqs  Rb/Mos  *  w^ere  Vvort  *s  me^sured  heat-transfer 
rate.  The  quantity  %vorf  -  q  was  determined  by  reducing  the  heat-transfer 
rate  measured  on  the  largest  model  of  the  test  by  the  square  root  of  the 
ratio  of  the  radii  and  by  the  ratio  of  the  heat-transfer  rate  for  the  largest 
model. to  the  heat-transfer  rate  to  this  model  for  the  high  Reynolds  number 
test. 

This  technique  is  used  to  eliminate  or  reduce  random  and  systematic 
errors  occurring  in  the  measurement  of  the  flow  conditions  and  in  the  model 
instrumentation.  The  use  of  this  technique  accounts  for  the  smaller  scatter 
of  the  PIBAL  data  relative  to  the  C.A.L.  data.  As  has  been  indicated,  the 
present  data  is  reduced  directly  from  the  measured  heat-transfer  rate  and 
the  measured  and  calculated  flow  conditions.  To  compare  the  two  sets  of 
data  with  Cheng's  theory  in  terms  of  absolute  values  of  heat-transfer  co¬ 
efficient,  a  zero- vorticity  heat-transfer  rate  has  been  calculated  for  the 

Z  6 

PIBAL*  test  conditions  using  the  theory  of  Fay  and  Riddell  .  The  details 
are  given  in  the  Appendix. 

Considering  the  degree  to  which  the  test  flow  conditions  differ  from 
the  assumptions  of  the  theory  (e.  g.  ,  Tw /T0  =  0.46  for  the  PIBAL  experi¬ 

ments  and  Tw/to  =  0  for  the  theory),  the  agreement  over  such  a  wide  range 
of  Reynolds  numbers  is  remarkably  good.  Furthermore,  the  theory  shows, 
Figure  11,  that  data  obtained  at  different  flow  conditions  can  be  correlated 
in  terms  of  parameters  such  as  A,  K*  or  €*Res  • 
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FUTURE  SHOCK  TUNNEL  DEVELOPMENT  AND 
LOW-DENSITY  RESEARCH 


The  present  experiments  have  demonstrated  the  application  of  the 
shock  tunnel  to  hypersonic,  rarefied- gasdynamic  research.  To  determine 
the  over-all  suitability  of  the  shock  tunnel  in  this  research  field,  a  thorough 
investigation  has  been  made  of  the  requirements  for  a  hypersonic,  losv-density 
wind  tunnel.  This  study  has  considered  instrumentation,  simulation  require¬ 
ments,  and  shock  tunnel  capabilities.  A  brief  review  of  these  items  will  be 
given  here. 

Instrumentation  Developments 

Quite  naturally  one  of  the  first  matters  of  concern  is  the  instrumen¬ 
tation,  since  it  is  of  little  value  to  produce  a  hypersonic  rarefied- gas  flow 
if  no  measurements  can  be  made.  The  present  heat-transfer  instrumentation 
will  sense  heat- transfer  rates  of  less  than  1  B.  T.  U.  / sq  ft/ sec.  This  is 
sufficient  to  measure  the  heat  transfer  about  one  foot  from  the  nose  of  a 
slender  cone  at  a  density  altitude  of  ZOO,  000  ft  at  Mach  5  and  325,  000  ft  at 
Mach  20.  Thus,  no  additional  development  of  heat-transfer  instrumentation 
is  needed  for  experiments  within  this  flight  envelope. 

The  pressure  transducers  now  in  use  have  sensitivities  sufficient  to 

-5 

measure  a  static  pressure  of  4  x  10  atm  which  corresponds  to  an  altitude 
of  330,000  ft.  These  transducers  are  modified  Altec-Lansing  microphones, 
Type  21BR- 180-2,  which  are  of  the  variable  capacitance  type.  In  addition, 
small- size  piezoelectric  transducers  having  comparable  sensitivity  and  a 

5 

much  lower  acceleration  sensitivity  are  currently  being  developed  at  CAL  . 

Present  force- sensing  instrumentation  is  highly  refined  as  are  the 
heat  transfer  and  pressure  techniques  and  is  adequate  for  measuring  lift, 
drag,  and  pitching  moment  on  blunt  bodies.  Further  improvements  can  be 
expected  to  provide  sufficient  sensitivity  to  extend  these  measurements  to 
slender  bodies.  The  free-flight  technique  reported  in  References  47  and  48 
can  also  be  used.  Thus,  instrumentation  techniques  are  available  for  making 
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heat  transfer,  pressure  and  force  measurements  in  hypersonic,  rarefied 
air  flow  in  a  shock  tunnel. 

Simulation  Requirements 

Next,  consideration  has  been  given  to  the  flow  quantities  and  parameters 
that  must  be  either  duplicated  or  simulated  in  a  low-density,  hypersonic  wind 
tunnel.  Among  the  more  obvious  of  these  are  Mach  number,  Reynolds  number, 
Knudsen  number,  shock-wave  density  ratio,  and  specific-heat  ratio.  Each 
experiment  requires  some  analysis  of  the  flow  phenomena  to  determine  the 
important  parameters  or  conditions  requiring  duplication.  In  general, 
however,  a  hypersonic,  low-density  wind  tunnel  should  be  capable  of  simul¬ 
taneously  achieving  both  high  Mach  numbers  and  hypervelocities  over  a  wide 
range  of  Reynolds  numbers.  Again,  the  shock  tunnel  has  a  proven  capability 
in  this  respect.  The  subject  of  simulation  requirements  and  shock  tunnel 
capabilities  is  covered  in  more  thorough  detail  in  References  5  and  49. 

Two  very  important  phenomena  that  may  be  encountered  in  hypersonic, 

low-density  wind  tunnels  are  thermodynamic  nonequilibrium  (or  relaxation 

46 

effects)  in  the  nozzle  flow  and  thick  nozzle- wall  boundary  layers.  Non¬ 
equilibrium  effects  can  be  avoided  by  restricting  the  stagnation  temperature 
to  values  below  which  vibration,  dissociation,  ionization  and  nitric-oxide 
formation  occur.  This,  however,  directly  limits  the  flow  velocity  to  values 
of  7000  ft/sec  or  less.  The  attainment  of  higher  speeds  requires  stagnation 
temperatures  of  2000PK  or  greater.  At  high  stagnation  temperatures,  non¬ 
equilibrium  effects  can  become  very  significant,  particularly  if  low  stagnation 
pressures  are  used.  It  appears  necessary,  therefore,  when  operating  at  high 

stagnation  temperatures,  also  to  operate  at  high  stagnation  pressures  to 

50 

minimize  nonequilibrium  effects  .  Low  Reynolds  numbers  are  then  obtained 
naturally  from  the  expansion  of  the  flow  to  very  high  Mach  numbers  and  low 
pressures  and  temperatures.  This  is  precisely  the  technique  that  must  be 
used  to  duplicate  flight  conditions  as  has  been  indicated  in  References  5  and 
49.  The  shock  tunnel  is  extremely  well- suited  for  this  type  of  operation, 
and  its  short  test  duration  is  a  distinct  advantage  for  high  stagnation  pressure 
and  temperature  operation. 
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V 


It  is  well  known  that  the  nozzle-wall  boundary-layer  thickens  as  the 

stagnation  pressure,  and  hence  free-stream  Reynolds  number,  is  decreased. 

51 

Boundary-layer  thickness  also  increases  strongly  with  Mach  number 
Consequently,  hypersonic  wind  tunnels  are  generally  so  plagued  by  such 
thick  boundary  layers  that  in  many  cases  only  a  small  core  of  good  flow 
exists  along  the  nozzle  centerline.  In  view  of  this  characteristic,  it  seems 
desirable  that  low-density,  hypersonic  nozzles  be  relatively  large,  say  on 
the  order  of  several  feet  or  greater  in  diameter.  Nozzles  of  this  size  would 
require  a  very  large  vacuum  pumping  system,  if  they  were  to  operate  con¬ 
tinuously.  Here,  the  shock  tunnel  possesses  a  favorable  feature  in  its 
intermittent  operation  because  extremely  high- capacity  pumping  systems 
are  not  required. 

Shock  Tunnel  Developments 

These  considerations  of  the  requirements  of  a  hypersonic,  low-density 
wind  tunnel  have  led  to  the  conclusion  that  the  shock  tunnel  possesses  many 
features  which  make  it  well- suited  for  rarefied- gasdynamic  research.  It 
should  be  mentioned,  however,  that  the  foregoing  discussion  is  primarily 
intended  to  bring  attention  to  the  advantages  of  the  shock  tunnel  and  is  not 
meant  to  imply  that  it  is  the  only  facility  suitable  for  research  in  hypersonic, 
rarefied  gasdynamics. 

In  view  of  the  apparent  advantages  possessed  by  the  shock  tunnel  as 
a  hypersonic,  low-density  facility,  CAL  has  undertaken  construction  of  a 
shock  tunnel  having  a  large-scale  nozzle  capable  of  expanding  the  air  from 
high  stagnation  pressures  and  temperatures  to  high  Mach  numbers  over  a 
wide  range  of  Reynolds  numbers.  This  tunnel  will  have  a  6-foot  diameter 
test  section  and  will  operate  at  stagnation  temperatures  up  to  8000°K  and 
stagnation  pressures  up  to  2000  atm.  These  conditions  will  permit  the 
attainment  of  hypervelocity  low- density  flow  with  minimized  nonequilibrium 
effects.  The  Mach  number  range  will  be  from  about  15  to  25,  and  the  test 
flow  duration  will  be  about  5  milliseconds.  Based  upon  the  success  attained 
during  the  present  experiments  utilizing  a  shock  tunnel  that  was  never  designed 
for  low-density  operation,  the  new  6-foot  shock  tunnel  should  prove  to  be  an 
extremely  useful  tool  for  research  in  hypersonic,  rarefied  gasdynamics. 
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SUMMARY 


The  experimental  results  presented  here  have  shown  that,  for  a  two- 
dimensional  stagnation  point  (transverse  cylinder),  the  boundary-layer 
theory  of  Fay  and  Riddell  predicts  the  heat  transfer  very  well  even  when  the 
flow  behind  the  shock  wave  is  fully  viscous.  No  increase  in  heat  transfer 
due  to  viscous-layer  effects  could  be  detected  on  the  transverse  cylinder. 

The  results  of  tests  of  an  axisymmetric  stagnation  point  (hemisphere 
cylinder)  indicate  that  vorticity-interaction  and  viscous-layer  effects  increase 
the  heat  transfer  significantly  above  that  predicted  by  boundary-layer  theory. 
A  new  theory  by  H.  K.  Cheng  at  CAL  has  been  used  to  correlate  the  hemi¬ 
sphere  data.  The  agreement  is  quite  good  over  the  entire  range  of  the  data 
which  extends  from  the  vorticity-interaction  regime  through  the  incipient 
merged  layer  regime. 

It  is  clear  that  additional  experiments  at  lower  Reynolds  numbers 
(higher  Knudsen  numbers)  are  needed  to  establish  the  transition  from  con¬ 
tinuum  to  free-molecule  theory.  The  new  6-foot  shock  tunnel  now  under 
construction  at  CAL  should  permit  extension  of  the  range  of  the  experiments 
to  the  free-molecule  flow  regime. 
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APPENDIX 


The  procedure  used  to  convert  the  PIBAL  data  from  the  form  of  pre¬ 
sentation  in  Reference  10  to  the  present  form,  Figure  11,  is  described  in 
this  Appendix.  The  approximations  and  assumptions  involved  are  indicated. 

The  experimental  data  are  presented  by  Ferri  et  al^  as  the  ratio 
Ivort  / *  vort  =  0  versus  Reynolds  number.  The  quantity  <f.vort  is  the 
measured  heat-transfer  rate.  The  quantity  £  vori  »  0  was  determined  experi¬ 
mentally  from  the  largest  diameter  model  for  which  the  vorticity  effects  are 
zero  and  the  square  root  of  the  ratio  of  the  radii  of  the  models.  The  Reynolds 
number  is  defined  as 


Re 


os 


03 


J^OS 


(Al) 


where  the  quantities  ,  HQS  and  /*os  are  the  density,  enthalpy  and  viscos¬ 
ity,  respectively,  evaluated  at  stagnation-point  conditions  behind  the  bow 
shock.  is  the  body  radius.. 

The  experimental  conditions  are  free- stream  Mach  number  of  8  and 
stagnation  temperatures  of  890°K  ( 1600#R)  and  1280°K  (2300°R).  The  range 
of  stagnation  pressures  is  not  specified;  however,  they  can  be  determined 
from  the  quoted  values  of  Reynolds  number.  The  ratio  of  wall-to- stagnation 
enthalpies  is  taken  as  0.  46  and  0.  34  for  the  890°K  and  1280°K  experiments, 
respectively. 

The  Reynolds  number  of  Reference  10,  R*0S>  is  related  to  the  Reynolds 
number  behind  the  shock  of  the  present  experiments  and  Cheng1  s  theory, 

,  as  follows 

*  A, 


-  €  /*0-H-/ho)  (fa)  (A2) 


The  notation  of  the  present  report  rather  than  Reference  10  is  used  here. 


31 


where  6  -/>*//>,>  HQ  -  H0f  and  ju s  has  been  taken  equal  to  which  is  a 

good  approximation  at  these  flow  conditions. 

To  express  the  PIBAL  data  given  as  %vorf  vort  •  oveTSUS 

=  —*RmJ?0S  in  terms  of  CH  =  9//>ooV~(H0-Hw)  and  4«  7 r  Pr/ - 
where  K2  =  g'  Ret  (• -Zi.  \  t  the  procedure  is  as  follows: 

V/"#  Tos  J 

(1)  is  calculated  from  Eq.  (A2) 

(2)  A  is  calculated  from  Eqs.  (Z)  and  (3) 

Z  3 

(3)  =  o  *s  ca^culate<i  using  the  theory  of  Fay  and  Riddell 

(4)  Cu  is  taken  as  ( ^°rt — \  C u  , 

H  \*vort-o)  "vort-O 

The  PIBAL  data  plotted  in  Figure  11  have  been  obtained  in  this  manner.  It 
has  been  necessary  to  assume  M =  8  and  T0  =  890°K  and  1Z80°K  exactly 
for  all  the  data  because  the  variation  of  M &  with  Reynolds  number  and  the 
repeatability  of  T0  for  different  tests  are  not  given. 
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FREE  MOLECULE 
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nnntr  ilv  u«g  i  uwn  «L«mi.v  ■  ■  ■ —  - .  _  .... 

COOLED  BLUNT  BODY  FLYING  AT  HYPERSONIC  SPEEDS  (ADAPTED  FROM  REF.  18) 


NCH  CAL  HYPERSONIC  SHOCK  TUNNEL 


CALCULATES  STAOUTIOR  PRESSURE,  fi°calc  ( ?si “) 

Figure  3  CORRELATION  OF  MEASURED  AND  CALCULATED  STAGNATION  PRESSURES 
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BASED  OB 


STAGNATION  PRESSURE  -  T>.  (psU) 


KNUDSEH  HO.  -K* 


Flgur*  7  HEMISPHERE  CYLINDER  HEAT  TRAHSFER  MODEL 


39 


Figure  8  STAGNATION-POINT  HEAT  TRANSFER  TO  A  TRANSVERSE  CYLINDER 

CORRELATED  WITH 
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HEAT  TRANSFER  COEFFICIENT 
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Figure  9  STAGNATION-POINT  HEAT  TRANSFTR  TO  A  TRANSVERSE  CYLINDER 


CORRELATED  WITHcffe, 
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FFEE  MOLECULE  FLOW 


Figure  10  STAGNATION  POINT  HEAT  TRANSFER  TO  A  HEMISPHERE  CYLINDER 
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